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We explain the effect of dark matter (flat rotation curve) using modified gravitational 
i dynamics. We investigate in this context a low energy limit of generalized general rel- 

J.^ ' ativity with a nonlinear Lagrangian C <x R n , where R is the (generalized) Ricci scalar 

and n is parameter estimated from SNIa data. We estimate parameter f3 in modified 
gravitational potential V(r) oc — — (1 + (ir~)^)- Then we compare value of f3 obtained 
from SNIa data with j3 parameter evaluated from the best fitted rotation curve. We 
find f3 ~ 0.7 which becomes in good agreement with an observation of spiral galaxies 
rotation curve. We also find preferred value of f2 mj o from the combined analysis of su- 
pernovae data and baryon oscillation peak. We argue that although amount of "dark 
energy" (of non-substantial origin) is consistent with SNIa data and flat curves of spiral 
galaxies are reproduces in the framework of modified Einstein's equation we still need 
substantial dark matter. For comparison predictions of the model with predictions of 
the ACDM concordance model we apply the Akaike and Bayesian information criteria 
of model selection. 

Keywords: Dark Energy; Dark Matter; Modified Gravity. 



1. Introduction 

Different astronomical observations [1,2] are pointing out that our Universe be- 
comes, at present time, in accelerating phase of expansion. In principle, there are 
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two quite different approaches to explain this observational fact. In the first ap- 
proach (which can be called substantial) it is assumed that universe is filled by 
mysterious perfect fluid violating the strong energy condition p x + 3p x > 0, where 
p x and p x are, respectively, the energy density and the pressure of this fluid. The 
nature as well as origin of this matter, called dark energy, is unknown until now. 
Among these approaches have appeared concordance ACDM model, which predicts 
that baryons contribute only about 4% of the critical energy density, non-baryonic 
cold dark matter (CDM) about 25% and the cosmological constant A (vacuum en- 
ergy) remaining 70%. Although ACDM model fits well SNIa data [3,4] this model 
offers only description of the observations not their explanation. From the method- 
ological view point the conception of mysterious dark energy seems to be effective 
physical theory only and motivates theorists for searching of alternative approaches 
in which nature of dark energy will be known at the very beginning. 

In the first approach it is assumed that Einstein's theory of general relativity 
is valid which reduces in practice (after assuming Robertson- Walker symmetry of 
space slices) to the case of Friedman - Robertson - Walker models. Nevertheless, 
theoretically it is not a'priori excluded the possibility of cosmology based on some 
extension of Einstein's general relativity. In this paper we consider such particular 
cases. 

On the other hand, there are alternative ideas of explanation, in which instead 
of dark energy some modifications of Fricdmann's equation are proposed at the very 
beginning. In these approaches some effects arising from new physics like brane cos- 
mologies, quantum effects, nonhomogencitics effects etc. can mimic dark energy by 
a modification of Fricdmann equation. Freese & Lewis [5] have shown that contri- 
butions of type p n to Fricdmann's equation ZH 2 = p e s, where p c s is the effective 
energy density and n is a constant, may describe such situations phenomenolog- 
ically. These models (called by their authors called the Cardassian models) give 
rise to acceleration, although the universe is flat, contains the usual matter and 
radiation without any dark energy components. This models have been tested by 
many authors (see for example [6,7,8,9,10,11,12]). What is still lacking is a funda- 
mental theory (like general relativity) from which these models can be derived after 
postulating Robertson Walker symmetry. 

In this paper we shall consider such particular type of generalization of Einstein's 
general relativity in which Lagrangian is proportional to R n , where R is generalized 
Ricci skalar. In particular, Einstein's general relativity is recovered if we put n = 
1. This theory is part of the larger class of so-called f(R) gravity, i.e. theories 
derived from gravitational Lagrangians that are analytical (usually polynomial) 
functions of R. (see e.g. [13,14,15,16,17,18,19,20,21]. In this approach (we called 
it non-substantial) instead of postulating mysterious dark energy it is assumed 
some extension of general relativity. Then effect of acceleration appeared naturally 
as a dynamical effect of the model. For modified gravity one can find Newtonian 
potential in non-relativistic limit and ask about possibility to explain flat rotation 
curve of spiral galaxies - major evidence for dark matter in the universe [22,23,24]. 
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(However, see also [25] for non flat rotational curves.) The main goal here is to 
explore power of this particular generalization of gravity in the context of dark 
energy and dark matter problems. We argue that although cosmology with modified 
Lagrangian C oc R n can explain " dark energy problem" but baryon oscillation test 
distinguishes value of density parameter of mater to be equal Q m ,o — 0.3, i.e. 
problem of dark matter is yet not solved within the framework of f(R) theories. We 
also demonstrate that models under considerations can reproduce rotation curves 
of spiral galaxies. 

The structure of the paper is as follows. In section II we define class of cosmolog- 
ical models of essential theory of gravity with lagrangian proportional to the Ricci 
scalar. Section III is devoted to analyze constraints on model parameters from SNIa, 
baryon oscillation peak and CMB shift. In section IV we investigate problem of ro- 
tation curves of spiral galaxies. Section V summarizes our results and formulates 
general conclusion that models modified gravity which are based on generalized 
lagrangians £ oc R n and Palatini formalism although solve the acceleration and flat 
rotation curves problems still favor fl m .o — 0.3. 

2. Cosmological models of nonlinear Palatini gravity 

Let us consider the simplest cosmological model of generalized Einstein's theory 
of gravity with Lagrangian C = f{R) which is function of the (generalized) Ricci 
scalar. The action is assumed in the form 

A = A grav + A mat = J (^U)f(R) + 2L mat («f))d i x (1) 

where f(R) oc R n and n is a constant. We also assumed that dynamical equation 
determining evolution of the cosmological model can be derived from the action 
through the Palatini formalism in which both metric g and symmetric connection 
F are regarded as an independent variables. Thus R = R(g, F) = g^R^iT) denotes 
generalized Ricci scalar (see e.g. [17] for details). 

Because of homogeneity and isotropy of the surface t = const is assumed, t - 
being a global cosmic time, we choose Robertson- Walker metric i.e. 

ds 2 = -dt 2 + a 2 (t) \ x ^ kr2 dr 2 + r 2 (d6 2 + sin 2 (0)^ 2 )] . (2) 

where k ~ 0, ±1 is curvature index, r, 9, (f> are usual spherical coordinates. Some 
properties of these theories (in the Palatini framework) have been already investi- 
gated by Capozziello et al. [26]. It has been demonstrated that under two popular 
choices f(R) oc R n and f(R) oc In R both models provide well fits to the SNIa data. 

Here we consider matter content in the form of perfect fluid which satisfy the 
conservation condition: 

^ = -3H(p + p) (3) 

where H = ^ (In a) is Hubble'a function. For convenience we assume simple form of 
equation of state (E.Q.S). After adopting Palatini formalism field equation reduces 
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to the ordinary second order differential equation which admits first integral in the 
form: 



H 2 = F(a). 



(4) 



The first integral (4) is usually called (generalized) Friedman equation. This is the 
first order differential equation in which right hand side is determined by matter 
content and curvature. Due to simple relation between the scale factor and the 
rcdshift z (1 + z = — ) formula (4) can be written in the form H 2 = F(z) (ao 
denotes the present value of the scale factor which corresponds to the redshift 
z = 0) In the system filled by both dust matter and radiation the function F(a) 
takes the following form: 



4n(2 - n)KT] rad 
3/3(n-3) 2 

where £ grav = — — 



H 2 



K Vdust 



2n 



3(3 - n) 



Kr ]dust 



+ 



2n 



(n-3) 



(5) 
0. If 



2 h^R n yfg, P is dimensional constant, w dus t = V dust I P dust 
n = (3 = 1 and r\ ra d = then the classical FRW dust filled model is recovered. 
Let us formulate some important remarks: 

(1) the formula (5) contains many nonphysical parameters which can be replaced by 
dimcnsionlcss density parameters fli defined for each additive contribution to 
the r.h.s. of (4). This in turn can be treated as a (fictitious or real) component of 
some effective energy density. Density parameter fij is defined for each energy 
component in standard way fli = pi/3HQ where Hq is present value of the 
Hubble function and pi is energy density of i — th fluid. 

(2) for our further analysis it is useful to separate this contribution on the r.h.s. of 
(4) which represents real dust matter scaling like a~ 3 from the non-substantial 
effects of generalized Lagrangian (related with n-parameter). Then our basic 
formula can be rewritten to the new more suitable form: 



H 



2n 



(3-n) 



-CI 



nonLO 



+n 1 , Q (i + z)^ n{2 



Tl) 3(l-n) 

-2-"nonJ,0(l + Z) ™ 



(n-3) 

where parameter fl nO ni,0 is determined from the constraint H(z = 0) = H 



(6) 



2n 



nonl,0 — 



(3-n) 



4n (2 



(n-3)' 



(7) 



Here k = is assumed for simplicity (for more general formulas see [27]). 



One can check that in the case of n = 1 one obtains Einstein de Sitter model 
filled with matter and radiation. The basic formula (6) will be suitable in the 
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next section for providing priors on fi m ,o which can be obtained from independent 
extragalactic measurements or baryon oscillation peak. 

From formula (6) on can drive a few conclusions. The first is that rejection of Q ri o 
in (6) doesn't eliminate automatically the dust term. On the other hand substitution 
Vdust — give rise to rejection of the second (radiation) term automatically. The 

3(1 — n) 

next observation arising from (6) is that term r2 nfm z.o(l + z ) ™ pl&ys a role of 
lapse function. Therefore one can re-scale original cosmological time following the 
rule: 1 1— > r: dr 2 — flnoni,o(l + z) s dt 2 and then obtain, after re-scaling density 
parameters, a flat model which is dynamically equivalent to the flat FRW model 
with : TO o = ^m.o^rrr an d £Vo = Slr.o i ?} 2 ~v) ■ Therefore, the exact solutions are 
well known in the form of t = t(Q TOi o, Q r ,o, z )- 

It is also worth to notice that equation (6) is equivalent to 



H ,2 



H, 



F{a) = or ±J-+V(a) = 0, (8) 



where ' = £p, dt\Ho\ = dr, V(a) = —p e ff\ — F ^ a . Due to particle like 
representation of the dynamics in the form (8) it is possible its investigation in 
terms of what H.-J. Schmidt calls classical mechanics with the lapse function [28]. 



3. Observational constraints on modified gravity parameters 

Within the framework of modified gravity, the acceleration originates from non- 
substantial contribution arising from curvature modification. This gives rise to neg- 
ative effective pressure and leads to self accelerating cosmology. 

3.1. Constraining model parameters from SNIa data 

The fundamental test for parameters of cosmological model is based on the lumi- 
nosity distance as a function of red-shift di(z) (the so-called Hubble diagram) 

^> = (1+s) i^K^f^))- o 

where S7fc o = —-m and T{x) = sinh(x), x, sin(x) for k = — 1,0, +1 respec- 
tively. For distant SNIa relation between luminosity distance <ii, absolute magni- 
tude M and directly observed their apparent magnitude m has the following form: 

H = m-M = 51og 10 d L + 25 = 51og 10 L> L + M , (2) 

where M. = — 51og 10 H a + 25 and D L = H dL. 

The goodness of fit is characterized by the parameter 

/..thcor ,,obs\2 
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where /i° hs is the measured value, /x* heor is the value calculated in the model under 
consideration, and Oi is the total measurement error. Assuming that supernovae 
measurements come with uncorrelated Gaussian errors, one can determine the like- 
lihood function C oc exp(— x 2 /2). The Probability Density Function (PDF) of cos- 
mological parameters [1] can be derived from Bayes' theorem. Therefore, one can 
estimate model parameters by using a minimization procedure. It is based on the 
likelihood function as well as on the best fit method minimizing \ 2 ■ 

In our analysis we used two samples of supernovae. One of them is "Gold" Riess 
et al. sample of 157 SNIa [3]. Second one is the sample of 115 supernovae compiled 
recently by Astier et al. [4]. This latest sample of 115 supernovae is our basic sample. 

For statistical analysis we have restricted the parameter Q m .o to the interval 
[0,1] and n to [—10.0,10.0] (except n — and additionally n = 3 for w = 0). 
Moreover, because of the singularity at n — 3, w — we have separated the cases 
n > 3 and n < 3 for w = in our analysis. Please note that fl nO ni,0 is obtained 
from the constraint H(z = 0) = H - The results of two fitting procedures performed 
on Riess and Astier samples with different prior assumptions for n are presented 
in Tables 1 and 2. In the Table 1 the values of model parameters obtained from 
the minimum the \ 2 are given, whereas in Table 2 the results from marginalized 
probability density functions are displayed. The best fit (minimum % 2 ) gives n ~ 2.6 
with the Astier et al. sample versus n ~ 2.1 with the Gold sample. In Figure 1 we 
present Probability distribution obtained with the Astier sample for the parameters 
Q m ,o an d n for non-linear gravity model, (case n < 3 marginalized over the rest 
of parameters). Please note that from Fig. 1 we obtain a very weak dependence of 
PDF on the matter density parameter if only f2 mj o > 0.05. Because bounce is a 
generic features of presented models for n > 2 [27] it is interesting to calculate from 
observational data probability that value of n paprameter is greated from two. We 
find that P(n > 2) ~ 0.99. It means that bounce is strongly favored over big-bang 
scenario like to in loop quantum gravity for example [32]. The Fig. 2 and Fig. 3 
shows likelihood contours on the plane (il m}0 ,n) obtained (from fits to the SNIa 
data and baryon oscilation peak test respecitvely) , obtained for non- linear gravity 
model, for the case n < 3 marginalized over M. 

Most popular are the Akaike information criteria (AIC) [33] and the Bayesian 
information criteria (BIC) [34]. We use this criteria to select model parameters 
providing the preferred fit to data. 

One of the important problem of modern observational cosmology is the so- 
called degeneracy problem: many models with dramatically different scenarios agree 
with the present day observational data. Information criteria for model selection 
[29] can be used, in some subclass of dark energy models, in order to overcome this 
degeneracy [30,31]. Most popular are the Akaike information criteria (AIC) [33] and 
the Bayesian information criteria (BIC) [34]. We use this criteria to select model 
parameters providing the preferred fit to data. 
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Table 1. The flat non-linear gravity model with dust and radiation. Results of statistical analysis 
performed with Astier et al. and Gold Riess et al. samples of SNIa as a minimum \ 2 best-fit. We 
separately analyzed the case n > 3 and n < 3 . 



sample 






n 


x 2 


Gold 


0.35 


< 0.01 


3.001 


180.7 


n < 3 


0.89 


0.23 


2.13 


181.5 


n > 3 


0.35 


< 0.01 


3.001 


180.7 


Astier 


0.01 


-1.47 


3.11 


108.7 


n < 3 


0.98 


0.08 


2.59 


108.9 


n > 3 


0.01 


-1.47 


3.11 


108.7 



Table 2. The flat non-linear gravity cosmological model with dust and radiation. The values of 
the parameters obtained from one dimensional PDFs calculated on the Astier et al. and the Gold 
Riess et al. SNIa samples. Because of the singularity at n = 3 we separately analyze the cases 
n > 3 and n < 3 . 



sample 




Qnonl,Q 


n 


Gold 


0.01 


0.26 


2.11 


n < 3 


1.00 


0.26 


2.11 


n > 3 


0.01 


-0.01 


3.001 


Astier 


0.01 


0.09 


2.56 


n < 3 


1.00 


0.09 


2.56 


n > 3 


0.01 


-0.01 


3.01 




0.00 0.20 e.40 0.60 0.80 i.ee i.5e 2.00 2.50 3.00 

fi m ,0 n 

Fig. 1. Probability distribution obtained with the Astier sample for the parameters f2 m o and n, 
marginalized over the rest of parameters. Non-linear gravity model, n < 3 . 



The AIC [33] is defined by 

AIC = -21n£ + 2d , 



(4) 
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1.80- 



1.60- 



O.OO 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00 



Fig. 2. The fiat non-linear gravity model with dust and radiation (n < 3). Likelihood contours on 
the (Q m fi,n) plane, marginalized over A4, obtained from fits to SNIa Astier et al. sample. 



where C is the maximum likelihood and d the number of model parameters. The 
best model, with a parameter set providing the preferred fit to the data, is that 
which minimizes the AIC. 

The BIC introduced by Schwarzc [34] is defined as 

BIC = -2\nC + d\nN , (5) 

where N is the number of data points used in the fit. While AIC tends to favor 
models with a large number of parameters, the BIC penalizes them more strongly, 
so the later provides a useful approximation to the full evidence in the case of no 
prior on the set of model parameters [35]. 

Please note that both values of information criteria have no absolute sense and 
only the relative values between different models are physically interesting. For the 
BIC a difference of 2 is treated as a positive evidence (6 as a strong evidence) 
against the model with larger value of BIC [36,37]. If we do not find any posi- 
tive evidence from information criteria, the models are treated as identical, while 
eventually additional parameters are treated as not significant. Therefore, the in- 
formation criteria offer a possibility to introduce a relation of weak ordering among 
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Table 3. Results of AIC and BIC performed on the Astier versus the Gold Riess samples of SNIa. 



sample 


AIC 


BIC 


ACDM Gold 


179.9 


186.0 


ACDM Astier 


111.8 


117.3 


Non-Lin. Grav. Gold 


186.6 


195.8 


Non-Lin. Grav. Astier 


114.7 


122.9 



considered models. 

In the Table 3 the value of AIC and BIC for the ACDM and the non-linear 
gravity models are presented. Note that for both samples we obtain with AIC and 
BIC for the ACDM model smaller values than for non-linear gravity. Most interst 
is using a Bayesian framework to compare the cosmological models, because they 
automatically penalize models with more parameters to fit the data. Based on these 
simple information criteria, we find that the SNIa data still favor the ACDM model, 
because under a similar quality of the fit for both models, the ACDM contains less 
parameters. 

3.2. CMB shift parameter 

For stringent and deeper constraint on model parameters we include in our analysis 
the so called (CMB) "shift parameter". This parameter is defined as: 

R^^^T(y(z lss )) (6) 

where R = 1.716±0.062 [38] and z lss = 1089 [39]. The i?-parameter determines the 
angular scale of the first acustic peak through the angular distance to last scattering 
and physical scale of the sound horizon. It is insensitive with respect perturbations 
and are suitable to constrain model parameter. The region allowed by the analysis 
of (CMB) "shift parameter" the plane (Q TOj o,n) for non-linear gravity model (for 
the case n < 3) is presented on the Fig 4 (lower panel). 

We obtain for non-linear gravity model the values of the model parameters 
^m,o = 0.67 and n — 1.03 as a best fit. Please note that this area is not allowed by 
SNIa data. 

3.3. Baryon oscillation peak 

Recently Fairbarn and Goobar [40] used baryon oscillation peak detected in 
the SDSS Luminosity Red Galaxies Survey [41] as a independent test of Dvali- 
Gabadadze-Porrati (DGP) brane model. They used constraint for: 
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so that E(z) = H(z)/H and z x = 0.35 yield A = 0.469 ± 0.017. The quoted 
uncertainty corresponds to one standard deviation, where a Gaussian probability 
distribution has been assumed. These constraints could be used for fitting cosmo- 
logical parameters (see also [4,40,42]). 

Fairbarn and Goobar [40] showed that the joint constraints for both SNIA data 
and the baryon oscillations peak ruled out flat DGP model at the 99% confidence 
level. Analogical analysis can be performed for our model. We obtain for non-linear 
gravity model the values of the model parameters fi TO ,o — 0.28, £l n onl,o = 0.33 and 
n = 2.53 as a best fit. On the Fig. 3 we show the region allowed by the baryon 
oscillation test on the plane (f2 m;0 ,n) for non-linear gravity model with dust and 
radiation (for the case n < 3). 



3.00- 



2.80- 



2.60- 



2.40- 



2.00- 



1.80- 



1.60- 




O.OO 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00 



Fig. 3. The fiat non-linear gravity model with dust and radiation (n < 3). Likelihood contours on 
the (U m fi,n) plane, marginalized over A4, obtained from baryon oscillation peak test. 
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3.4. Combined SNIa, CMB shift and baryon oscillation 
constraints 

Now it is possible obtained constraints from both SNIa data and CMB shift and 
baryon oscillation peak. The results of our combined analysis are presented in the 
Fig. 4 (upper panel) On can see that the combination of three independent obser- 
vational constraints distinguish the value fl m ,o — 0.3 like for ACDM concordance 
model in which there is present substantial conception of dark matter. However 
please note that area allowed by CMB shift is excluded by area allowed by com- 
bined SNIa data and baryon oscillation peak because different value of n obtained 
in both cases. 



4. Flat rotation curves from cosmology in C oc R n theories 

Let us consider low energy limit of modified gravity with lagrangian L oc R n . 
For this aims it is useful to consider point like m in Schwarzchild - like metric 
(spherically symmetric). Then modified gravitational potential which corrected the 
ordinary Newtonian potential is of the form: 

T/(r)cx-i(l + (f)' 3 ) (1) 

where r c is characteristic parameter which crucially depends on the mass of the 
system and (3 = (3{n) [43]. 

Hence we can evaluate the rotation curve in the Newtonian limit of modified 
gravity mv ^ r ) = _ j n ^ ne previous section we estimate value of n ~ 2.6. Then 
we calculate (3 parameter from the formula: 

_ 12n 2 - In - 1 - V36n 4 + 12n 3 - 83n 2 + 50n + 1 
1 6n 2 + An - 2 ( ' 

obtained by Capozzielo et al. [43]. 

We obtain (3 ~ 0.7 which is close to (3 estimated for NGC 5023 ((3 = 0.714). 
Therefore we obtain that considered theory reproduce flat rotation curves of spiral 
galaxies. Moreover, the value of (3 parameter required to explain acceleration ex- 
pansion of the Universe give rise to correct peculiarities of observed rotation curve. 
Nevertheless note, that from investigation presented in previous section density pa- 
rameter for matter is close to £l m ,o = 0.3 rather than to value Q m ,o = 0.05 as can be 
expected if both effects of dark energy and dark matter has non-substantial nature 
i.e. (arises from modified gravity only). 

5. Conclusion 

In this paper we consider the simplest choice of f(R) theories with f(R) oc R n . The 
basic motivation is searching for fundamental theory of gravity capable to explain 
both dark energy and dark matter problems without referring to mysterious dark 
energy conception. For this aim we consider cosmology based on such a theory of 
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gravity and then we use different observational constraints on independent model 
parameters. We consider simple flat FRW model. It is intcgrablc in exact form after 
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re-parametrization of time variable. From estimation based on SNIa and BOP we 
obtain n > 2 which means that bouncing phase instead of big-bang singularity is 
generic features of such models. Because n > 2 new r parameter is monotonous 
function of cosmic time and acceleration epoch is transitional only phenomenon. In 
the future the universe decelerate which distinguish our model from ACDM one. 
Note that because for small value of scale factor a curvature effects are negligible in 
the comparison to other matter contribution, therefore, in the generic case big-bang 
singularity is replaced by bounce. 

Analysis of SNIa Astier data shows that values of \ 2 statistic are comparable for 
both ACDM and best fitted non-linear gravity model. For deeper analysis we use 
Akaike and Bayesian information criteria of model comparison and selection. We 
find these criteria still to favor the ACDM model over non-linear gravity, because 
(under the similar quality of the fit for both models) the ACDM model contains 
one parameter less. 

Moreover, we find that the effect of dark matter can be kinematically explained 
as a effect of nonlinear gravity with Lagrangian C oc R n . Parameter /3 required for 
explaining accelerated expansion of the universe give rise to correct peculiarities 
of observed rotation curve. However from baryon oscillation peak prior we still 
obtain Cl m ,o — 0.3 (instead of fi m ,o ~ 0.05 as we expected). Moreover, we find 
a disagreement between results obtained from CMB shift parameter analysis and 
that from joint SNIa and baryon oscillation peak. Finally, the substantial form of 
dark matter is still required. 
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